Self-organized transition to coherent activity in disordered media 



(N 

o 

(N 
G 

in 



O 
d 



> 



O 

(N 



X 



Rajeev Singh 1 , Jinshan Xu 2,3 , Nicolas Garnier 2 , Alain Pumir 2 and Sitabhra Sinha 1 
The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai 600113, India. 
2 Laboratoire de Physique, ENS de Lyon and CNRS, 46 Allee d'ltalie, 69007, Lyon, France. 
3 Department of Physics, East China Normal University, Shanghai, 20062, China. 

(Dated: January 6, 2012) 

Synchronized oscillations are of critical functional importance in many biological systems. We 
show that such oscillations can arise without centralized coordination in a disordered system of 
electrically coupled excitable and passive cells. Increasing the coupling strength results in waves 
that lead to coherent periodic activity, exhibiting cluster, local and global synchronization under 
different conditions. Our results may explain the self-organized transition in a pregnant uterus from 
transient, localized activity initially to system-wide coherent excitations just before delivery. 
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Rhythmic behavior is central to the normal function- 
ing of many biological processes [l[ and the periods of 
such oscillators span a wide range of time scales control- 
ling almost every aspect of life [2|-|5| . Synchronization of 
spatially distributed oscillators is of crucial importance 
for many biological systems [f|. For example, disrup- 
tion of coherent collective activity in the heart can result 
in life-threatening arrhythmia [7j. In several cases, the 
rhythmic behavior of the entire system is centrally orga- 
nized by a specialized group of oscillators (often referred 
to as pacemakers) Q as in the heart, where this func- 
tion is performed in the sino-atrial node However, 
no such special coordinating agency has been identified 
for many biological processes. A promising mechanism 
for the self-organized emergence of coherence is through 
coupling among neighboring elements. Indeed, local in- 
teractions can lead to order without an organizing center 
in a broad class of complex systems ■ 

The present work is inspired by studies of the pregnant 
uterus whose principal function is critically dependent on 
coherent rhythmic contractions that, unlike the heart, do 
not appear to be centrally coordinated from a localized 
group of pacemaker cells [ll|. In fact, the uterus re- 
mains quiescent almost throughout pregnancy until at 
the very late stage when large sustained periodic activity 
is observed immediately preceding the expulsion of the 
fetus [12}. In the USA, in more than 10 % of all preg- 
nancies, rhythmic contractions are initiated significantly 
earlier, causing preterm births fl3j . which are responsible 
for more than a third of all infant deaths [3] ■ The causes 
of premature rhythmic activity are not well understood 
and at present there is no effective treatment for preterm 
labor [l2|. 

In this paper we have investigated the emergence of 
coherence using a modeling approach that stresses the 
role of coupling in a system of heterogeneous entities. 
Importantly, recent studies have not revealed the pres- 
ence of pacemaker cells in the uterus [Hj]. The uterine 
tissue has a heterogeneous composition, comprising elec- 
trically excitable smooth muscle cells (uterine myocytes), 



as well as electrically passive cells (fibroblasts and inter- 
stitial Cajal-like cells [ICLCs]) [l6, 17 1. Cells are coupled 
in tissue by gap junctions that serve as electrical conduc- 
tors. In the uterine tissue, the gap junctional couplings 
have been seen to markedly increase during late preg- 
nancy and labor, both in terms of the number of such 
junctions and their conductances (by an order of magni- 
tude [HI), which is the most striking of all electrophys- 
iological changes the cells undergo during this period. 
The observation that isolated uterine cells do not spon- 
taneously oscillate [l5|, whereas the organ rhythmically 
contracts when the number of gap junctions increases, 
strongly suggests a prominent role of the coupling. The 
above observations have motivated our model for the on- 
set of spontaneous oscillatory activity and its synchro- 
nization through increased coupling in a mixed popula- 
tion of excitable and passive elements. While it has been 
shown earlier that an excitable cell connected to passive 
cells can oscillate [l9| , we demonstrate that coupling such 
oscillators with different frequencies (because of varying 
numbers of passive cells) can result in the system having a 
frequency higher than its constituent elements. We have 
also performed a systematic characterization for the first 
time of the dynamical transitions occurring in the het- 
erogeneous medium comprising active and passive cells 
as the coupling is increased, revealing a rich variety of 
synchronized activity in the absence of any pacemaker. 
Finally, we show that the system has multiple coexist- 
ing attractors characterized by distinct mean oscillation 
periods, with the nature of variation of the frequency 
with coupling depending on the choice of initial state as 
the coupling strength is varied. Our results provide a 
physical understanding of the transition from transient 
excitations to sustained rhythmic activity through phys- 
iological changes such as increased gap junction expres- 
sion [H. 

The dynamics of excitable myocytes can be described 
by a model having the form C m V e = —Iion{V e ,gi) where 
V e (mV) is the potential difference across a cellular mem- 
brane, C m (= 1 /iF cm~ 2 ) is the membrane capacitance, 
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Ii on (a*A cm -2 ) is the total current density through 
ion channels on the cellular membrane and Qi are the 
gating variables, describing the different ion channels. 
The specific functional form for I ion varies in different 
models. To investigate the actual biological system we 
have first considered a detailed, realistic description of 
the uterine myocyte given by Tong et al. [2l[. How- 
ever, during the systematic dynamical characterization 
of the spatially extended system, for ease of computation 
we have used the phenomenological FitzHugh-Nagumo 
(FHN) system Q which exhibits behavior qualitatively 
similar to the uterine myocyte model in the excitable 
regime. In the FHN model, the ionic current is given by 
lion = Fe(V e ,g) = AV e (V e - a)(l - V e ) - g, where g is 
an effective membrane conductance evolving with time 
as g = e(V e — g), a(= 0.2) is the excitation threshold, 
A(= 3) specifies the fast activation kinetics and e(= 0.08) 
characterizes the recovery rate of the medium (the pa- 
rameter values are chosen such that the system is in the 
excitable regime and small variations do not affect the re- 
sults qualitatively). The state of the electrically passive 
cell is described by the time-evolution of the single vari- 
able V p HI: V p = Fp{V p ) = K(V R - V p ), where the rest- 
ing state for the cell, V R is set to 1.5 and K(= 0.25) char- 
acterizes the time-scale over which perturbations away 
from V R decay back to it. We model the interaction 
between a myocyte and one or more passive cells by: 

V e = F e (V e ,g) + n p C r (V p - V e ), (la) 
V p = Fp(Vp)-C r (V p -V e ), (lb) 

where n p (= 1,2,.. .) passive elements are coupled to an 
excitable element via the activation variable V e , p with 
strength C r . Here, we have assumed for simplicity that 
all passive cells are identical having the same parame- 
ters V p R and K, as well as, starting from the same initial 
state. We observe that the coupled system comprising a 
realistic model of uterine myocyte and one or more pas- 
sive cells exhibits oscillations (Fig. [1] (a)) qualitatively 
similar to the generic FHN model (Fig. [TJ (b)), although 
the individual elements are incapable of spontaneous pe- 
riodic activity in both cases. In Fig. [TJ (a-b), the range 
of n p and excitable-passive cell couplings for which limit 
cycle oscillations of the coupled system are observed is in- 
dicated with a pseudocolor representation of the period 
(t) . We also look at how a system obtained by diffusively 
coupling two such "oscillators" with distinct frequencies 
(by virtue of having different n p ) behaves upon increas- 
ing the coupling constant D between V e (Fig. Q] (c)). A 
surprising result here is that the combined system may 
oscillate faster than the individual oscillators comprising 
it. 

To investigate the onset of spatial organization of pe- 
riodic activity in the system we have considered a 2- 
dimensional medium of locally coupled excitable cells, 
where each excitable cell is connected to n p passive cells 
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FIG. 1: (Color online) Oscillations through interaction be- 
tween excitable and passive elements. A single excitable el- 
ement described by (a) a detailed ionic model of an uterine 
myocyte and (b) a generic FHN model, coupled to n p pas- 
sive elements exhibits oscillatory activity (inset) with period 
r for a specific range of gap junctional conductances G m in 
(a) and coupling strengths CV in (b). The triangles (upright 
and inverted) enclosing the region of periodic activity in (b) 
are obtained analytically by linear stability analysis of the 
fixed point solution of Eq. (|lap . (c) Frequency of oscillation 
for a system of two "oscillators" A and B (each comprising an 
excitable cell and n p passive cells with = 1 and n p = 2) 
coupled with strength D. Curves corresponding to different 
values of CV show that the system synchronizes on increas- 
ing D, having a frequency that can be higher than either of 
the component oscillators, (d) Uterine tissue model as a 2- 
dimensional square lattice, every site occupied by an excitable 
cell coupled to a variable number of passive cells. 



[Fig. [TJ (d)], n p having a Poisson distribution with mean 
/. Thus, / is a measure of the density of passive cells rel- 
ative to the myocytes. Our results reported here are for 
/ = 0.7; we have verified for various values of / > 0.5 
that qualitatively similar behavior is seen. The dynamics 
of the resulting medium is described by: 



8Ve 

dt 



F e (V ei g) + n p C r {V p - V e ) + DV 2 V e , 



where D represents the strength of coupling between ex- 
citable elements (passive cells are not coupled to each 
other). Note that, in the limit of large D the behavior 
of the spatially extended medium can be reduced by a 
mean-field approximation to a single excitable element 
coupled to / passive cells. As / can be non-integer, n p 
in the mean-field limit can take fractional values [as in 
Fig. Q] (a-b)]. 

We discretize the system on a square spatial grid of 
size L x L with the lattice spacing set to 1. For most re- 
sults reported here L = 64, although we have used L up 
to 1024 to verify that the qualitative nature of the tran- 
sition to global synchronization with increasing coupling 
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FIG. 2: (Color online) Emergence of synchronization with in- 
creased coupling. Snapshots (first row) of the activity V e in a 
two-dimensional simulation domain (/ = 0.7, C r = 1, L = 64) 
for increasing values of coupling D (with a given distribution 
of n p ). The corresponding time-averaged spatial correlation 
functions C(r) are shown in the middle row. The size of the 
region around r = (at center) where C(r) is high provides 
a measure of the correlation length scale which is seen to 
increase with D. The last row shows pseudocolor plots indi- 
cating the frequencies of individual oscillators in the medium 
(white corresponding to absence of oscillation). Increasing D 
results in decreasing the number of clusters with distinct os- 
cillation frequencies, eventually leading to global synchroniza- 
tion characterized by spatially coherent, wavelike excitation 
patterns where all elements in the domain oscillate with same 
frequency. 



is independent of system size. The dynamical equations 
are solved using a fourth-order Runge Kutta scheme with 
time-step dt < 0.1 and a standard 5-point stencil for the 
spatial coupling between the excitable elements. We have 
used periodic boundary conditions in the results reported 
here and verified that no-flux boundary conditions do not 
produce qualitatively different phenomena. Frequencies 
of individual elements are calculated using FFT of time- 
series for a duration 2 15 time units. The behavior of 
the model for a specific set of values of /, C r and D is 
analyzed over many (~ 100) realizations of the n p distri- 
bution with random initial conditions. 

Fig- HI (first row) shows spatial activity in the sys- 
tem at different values of D after long durations (~ 2 15 
time units) starting from random initial conditions. As 
the coupling D between the excitable elements is in- 
creased, we observe a transition from highly localized, 
asynchronous excitations to spatially organized coherent 
activity that manifests as propagating waves. Similar 
traveling waves of excitation have indeed been experi- 
mentally observed in vitro in myometrial tissue from the 
pregnant uterus (23[. The different dynamical regimes 
observed during the transition are accompanied by an in- 
crease in spatial correlation length scale (Fig. [21 middle 
row) and can be characterized by the spatial variation of 
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FIG. 3: (Color online) (a) Different dynamical regimes of the 
uterine tissue model (for / = 0.7) in D — C' r parameter plane 
indicating the regions having (i) complete absence of oscilla- 
tion (NO), (ii) cluster synchronization (CS), (iii) local syn- 
chronization (LS), (iv) global synchronization (GS) and (v) 
coherence (COH). (b-c) Variation of (b) width of frequency 
distribution (ov) and (c) fraction of oscillating cells (/ OJC ) 
with coupling strength D for C T = 1 [i.e., along the broken 
line shown in (a)]. The regimes in (a) are distinguished by 
thresholds applied on order parameters (<7„), (fosc) and (F), 
viz., NO: (f 03C ) < 10~ 3 , CS: (a v ) > 10" 4 , LS: (er„) < 10" 4 
and (/ osc ) < 0.99; GS: </ osc ) > 0.99 and COH: (F) > 0.995. 
Results shown are averaged over many realizations. 



frequencies of the constituent elements (Fig. [21 last row). 
For low coupling (D = 0.1), multiple clusters each with a 
distinct oscillation frequency v coexist in the medium. As 
all elements belonging to one cluster have the same pe- 
riod, we refer to this behavior as cluster synchronization 
(CS). Note that there are also quiescent regions of non- 
oscillating elements indicated in white. With increased 
coupling the clusters merge, reducing the variance of the 
distribution of oscillation frequencies eventually result- 
ing in a single frequency for all oscillating elements (seen 
for D = 0.3). As there are still a few local regions of 
inactivity, we term this behavior as local synchronization 
(LS). Further increasing D (=0.4), a single wave traverses 
the entire system resulting in global synchronization (GS) 
characterized by all elements in the medium oscillating 
at the same frequency. Our results thus help in causally 
connecting two well-known observations about electrical 
activity in the pregnant uterus: (a) there is a remarkable 
increase in cellular coupling through gap junctions close 
to onset of labor [l8| and (b) excitations are initially in- 
frequent and irregular, but gradually become sustained 
and coherent towards the end of labor [ll[ . 

The above observations motivate the following order 
parameters that allow us to quantitatively segregate the 
different synchronization regimes in the space of model 
parameters [Fig. [3] (a)]. The CS state is characterized by 
a finite width of the frequency distribution as measured 
by the standard deviation, a v , and the fraction of oscil- 
lating elements in the medium, < f osc < 1. Both LS 
and GS states have o~ v — > 0, but differ in terms of f osc 
(< 1 in LS, ~ 1 in GS). Fig. [3] (b-c) shows the variation 
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FIG. 4: (Color online) (a) Variation of mean oscillation fre- 
quency v with coupling strength D in the uterine tissue model 
(/ = 0.7) for 400 different initial conditions at CV = 1. Con- 
tinuous curves correspond to gradually increasing D starting 
from a random initial state at low D, while broken curves 
(overlapping) correspond to random initial conditions chosen 
at each value of D. (b-c) Snapshots of activity in the medium 
at D = 1.5 for a random initial condition seen at intervals of 
ST = 5 time units, (d) Variation of the cumulative fractional 
volumes (f> of the basins for different attractors corresponding 
to activation patterns shown in (b-c) and (e), as a function of 
the coupling strength D. (e) Snapshots of topologically dis- 
tinct patterns of activity corresponding to the five attractors 
at D = 1.5 (broken line in (d)) when D is increased. 



of the two order parameters (ov) and (/ OS c) with the cou- 
pling D, ( ) indicating ensemble average over many real- 
izations. Varying the excitable cell-passive cell coupling 
C r together with D allows us to explore the rich variety 
of spatio-temporal behavior that the system is capable 
of [Fig. |3] (a)]. In addition to the different synchronized 
states (CS, LS and GS), we also observe a region where 
there is no oscillation (NO) characterized by f osc — > 0, 
and a state where all elements oscillate with the same 
frequency and phase which we term coherence (COH). 
COH is identified by the condition that the order param- 
eter F = max t [/ act (t)] — > 1 where f ac t(t) is the fraction 
of elements that are active (V e > a) at time t. In prac- 
tice, the different states are characterized by thresholds 
whose specific values do not affect the qualitative nature 
of the results. 

To further characterize the state of the system, we 
determined the mean frequency v by averaging over all 
oscillating cells for any given realization of the system. 
Fig. [4] (a) reveals that several values of the mean fre- 
quency are possible at a given coupling strength. When 
the initial conditions are chosen randomly for each value 
of the coupling (broken curve in Fig.@](a)), the mean fre- 
quency decreases with increasing D. On the other hand, 
v is observed to increase with D when the system is al- 



lowed to evolve starting from a random initial state at 
low D, and then adiabatically increasing the value of D. 
The abrupt jumps correspond to drastic changes in the 
size of the basin of an attractor at certain values of the 
coupling strength, which can be investigated in detail 
in future studies. This suggests a multistable attractor 
landscape of the system dynamics, with the basins of the 
multiple attractors shown in Fig. 2] (d) [each correspond- 
ing to a characteristic spatiotemporal pattern of activity 
shown in Fig. U (e)] having differing sizes. They repre- 
sent one or more plane waves propagating in the medium 
and are quite distinct from the disordered patterns of 
spreading activity (Fig. U (b-c)) seen when random ini- 
tial conditions are used at each value of D. We note that 
the period of recurrent activity in the uterus decreases 
with time as it comes closer to term (2pj in conjunction 
with the increase in number of gap junctions. This is 
consistent with our result in Fig. |4ja) when considering 
a gradual increase of the coupling D. 

Our results explain several important features known 
about the emergence of contractions in uterine tissue. 
Previous experimental results have demonstrated that 
the coupling between cells in the myometrium increases 
with progress of pregnancy [HI- This suggests that the 
changes in the system with time amounts to simulta- 
neous increase of D and C r , eventually leading to syn- 
chronization as shown in Fig. [3] (a) . Such a scenario is 
supported by experimental evidence that disruption of 
gap-junctional communication is associated with acute 
inhibition of spontaneous uterine contractions p4l |. The 
mechanism of synchronization discussed here is based on 
a very generic model, suggesting that our results apply 
to a broad class of systems comprising coupled excitable 
and passive cells (25l [26j . A possible extension will be to 
investigate the effect of long-range connections 27 1 . 

To conclude, we have shown that coherent periodic ac- 
tivity can emerge in a system of heterogeneous cells in a 
self-organized manner and does not require the presence 
of a centralized coordinating group of pacemaker cells. A 
rich variety of collective behavior is observed in the sys- 
tem under different conditions; in particular, for inter- 
mediate cellular coupling, groups of cells spontaneously 
form clusters that oscillate at different frequencies. With 
increased coupling, clusters merge and eventually give 
rise to a globally synchronized state marked by the gen- 
esis of propagating waves of excitation in the medium. 
Our model predicts that a similar set of changes occur in 
the uterus during late stages of pregnancy. 
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